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The thermal conductivity of pionic medium has been evaluated with the help of its 
standard expression from the relaxation time approximation, where inverse of pion re¬ 
laxation time or pion thermal width has been obtained from the imaginary part of pion 
self-energy. In the real-time formalism of thermal field theory, the finite temperature 
calculations of pion self-energy for na and tt p loops have been done. The numerical 
value of our thermal conductivity increases with temperature very softly, though at par¬ 
ticular temperature, our estimation has to consider a large band of phenomenological 
uncertainty. 
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1. Introduction 

In order to describe the expanding fire ball, produced at RHIC or LHC, the inclusion 
of dissipative effects in hydrodynamical and transport simulations has attracted 
remarkable attention in the recent years. Only the shear viscosity of the expanding 
matter is taken as a relevant dissipative coefficient by most of the communities, 
doing such kind of simulations. To explain the elliptic flow parameter, V 2 , extracted 
from data collected at RHIC and LHC, their investigations suggest that the matter 
is likely to have a very small ratio of shear viscosity to entropy density, rj/s. As a 
next attempt, the inclusion of bulk viscosity coefficient (£) in those simulations has 
been started in a few studies but the effect of thermal conductivity (k) on those time 
evolution pictures of expanding matter is largely ignored except the recent Refsp® 
trying to develop a model of the thermal conductivity. This is because the thermal 
conduction does not take place very well in the system, where net baryon density 
is approximately zero.® ° Most of the hydrodynamical simulations are focused on 
describing the baryon-free matter, which is expected to be created in the central 
rapidity region at the RHIC and LHC. However, thermal conductivity definitely 
demands its attention to describe the expanding fire ball at forward rapidities or at 
smaller collision energies, such as those at FAIR and in the low-energy runs at RHIC. 
Probably for same reason, the microscopic calculation of thermal conductivity for 
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Table 1. Numerical values of thermal conductivity k of some 
earlier works in the hadronic temperature domain T Ri 0.12 — 
0.17 GeV or near the transition temperature. 


Earlier works (adopted model or method) 

k (GeV 2 ) 

Ref.ES (NJL model) 

11.6-7.78 

RefES (DPQM model) 

4.5 

Ref® (Liquid-instanton model) 

0.12 

Ref! 11 ! (unitarization method) 

0.07-0.1 

Ref.® (CEA method) 

0.032-0.049 

Ref® (RTA method) 

0.016-0.02 

Ref.® (BAMPS, a Ri43 mb) 

2.59/cr & 0.023 

Ref® (CEA method) 

0.015-0.024 

Ref® (RTA method) 

0.016-0.024 

Ref® (CEA method) 

0.016-0.024 

RefPl (unitarization method) 

0.007-0.008 


the strongly interacting matter is not extensively explored to know its explicit 
temperature dependence as done rigorously for the other dissipative quantities viz. 
shear and bulk viscosities. To date, only a few studies have addressed calculation 
of thermal conductivity of the strongly interacting matter! 3 16 Their estimations 
collectively exhibit a band of predictions as shown in the Table [l| 

Here our interest is to calculate the thermal conductivity of pionic medium, 
which may also appear in between the chemical and kinetic freeze-out temperatures 
( T c h ~ 0.17 GeV to Tf « 0.12 GeV) at zero baryon chemical potential of the 
expanding fire ball during its late stage. Though the pionic constituent particles 
do not carry any baryon number but for the medium with finite pion chemical 
potential /i 7r , where total number of pion remain constant, may have a non-zero 
value of thermal conductivityP® We will make our estimation for = 0, which 
will help to compare with other relevant results. 


2. Formalism 

Let us start from the standard expression of thermal conductivity for pionic medium 
in the relaxation time approximation (RTA) 3 : 

(4 - 4) V(4){ 1 + M4)l . (!) 

where nfc(w^) = 1 /{e^ k — 1} is the Bose-Einstein (BE) distribution function of 
pion with = ( k 2 + m 2 ) 1//2 , and T^^ l/r ff , tv is pion relaxation time) is the 
thermal width of 7r mesons in the medium. Here /„• = 3 is pion isospin factor and 
hn is the exact heat function per particle for the pionic medium, which can be 
estimated from the ideal gas expression, 


_ /4G f °° dkk 4 

6tt 2 J 0 w£ 2 R 


/i-TT — (&-K T Ptt)/P n — Ts-tt/ Pit , 


(2) 
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Fig. 1. Pion self-energy diagram for ircr and np loops. 


where 

f d 3 k ( k 2 \ 

(3) 

and 

f d 3 k . . . 

Pn In J ^27r) 3 ' (^) 

are entropy density and number density respectively of the pionic medium at tem¬ 
perature T = 1/p. 

Now, the thermal width, for pion can be determined from the imaginary part 
of its one-loop self-energjt 1720 at finite temperature, by the relation: 

r 7 r (fc,r) = T ) 

M 

= - lmU n(nM) ( fc 0 = <4, fc, T)/m n ■ (5) 

M 


In the above equation, 11^^ is the retarded component of pion self-energy for 
7 tM loops, where M stands for a and p resonances. Feynman diagram of pion self 
energy for irM loops is shown in the Fig. ©, where k is the momentum of external 
7T line while l and u = k — l are momenta of internal 7r and M lines respectively. 

Following the effective hadronic model as adopted in earlier Refsf^^Ql^Il the a 
and p resonances in 7T7 t scattering have traditionally been introduced by using the 
effective Lagrangian densities for 7T7rer and nnp interactions: 


C = g P Pp ■ n x <9 m tt 


-m a TT ■ 7r <t, 


( 6 ) 


where the coupling constants (g a = 5.82 and g p = 6) are fixed from their experi¬ 
mental decay widths of a and p mesons in their 7T7 t channels.^ The negative sign 
in the second term for ann coupling is important to determine the relative phased 
With the help of the above Lagrangian densities ©, pion self-energy for irM loops 
(shown in Fig. [T]) with M = a and p can be derived in the real-time formalism of 
thermal field theory.^ As the pion pole (ko = k) is situated in the Landau cut 
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Fig. 2. Invariant mass distribution function, r 7 T( 7 T M){Mk) = Imll ^ nM y{M k )/m n for 7T(j (dotted 

line), 7r p (dashed line) loops and their total (solid line) at fixed values of k and T. Upper and 
lower panels consist the results with and without folding effect respectively. 


of the self-energy function, k) so the corresponding thermal width can 

simply be expressed as 17 M22 


r 7r(7rM)(fc,r) = Imn^ (wM) (fc 0 = w l,k,T)/m n 


1 


d 3 l 


L(lo 


-uf,l,k){ni{u i) 


-n u {u™)}8{k 0 + ui ~^u)\ ko=l 


(7) 


where ni and n u are BE distribution functions of 7r and M mesons respectively. 
Using the effective Lagrangian densities ©, one can obtain the vertex factors: 


L(k,l) = for M = a , 

= [k 2 (fc 2 — m 2 ) + l 2 (l 2 — to 2 ) — 2 {(k ■ l ) to 2 + k 2 Z 2 }], for M = p . 

trip 

( 8 ) 

A hadronic form factor F(l) = A 2 /(JY + A 2 ) with A = 1 GeV has been multiplied 
with each of the effective coupling constants, g a and g p to consider finite size effect 
of hadronic vertices. 


3. Results and Discussion 

Let us start our discussion of numerical results from the pion thermal widths, which 
inversely control the numerical strength of thermal conductivity of pionic medium. 
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Fig. 3. At fixed pion momentum (k = 0.5 GeV), the temperature dependence of (dotted 

line), (dashed line) and their total (solid line) are shown in the upper panel whereas lower 

panel demonstrates the corresponding mean free path contributions of two individual loops and 
their total. 



k (GeV) 


Fig. 4. Similar to Fig. for a fixed value of T against k- axis. 


In Fig. (0, the off-mass shell behavior of r w ( 7rM )(Mfe, k, T) is plotted against the 
invariant mass axis Mf,. before going to analyze the on-shell function M){M k = 
win ,k,T). The numerical values of T^^m) (Mk ,k,T) can be generated if we replace 
k 0 = {k 2 + m 2 } 1 / 2 by k 0 = {k 2 + M 2 } 1 / 2 in Eq. 0 for fixed values of pion 
momentum (k = 0.5 GeV) and temperature of the medium (T = 0.13 GeV). In 
the lower panel of Fig. 0, the dotted line is representing the Landau and unitary 
cut contributions for na loop in the regions (M*. = 0 to m CT — m n = 0.25 GeV) 
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and (Mfc = m a + m^ = 0.53 GeV to oo) respectively. Similarly, the dashed line is 
displaying the corresponding Landau and unitary cut contributions for np loop in 
the regions (M*, = 0 to m p — m^ = 0.63 GeV) and (.M*, = m p + m n = 0.91 GeV to 
oo) respectively. Owing to the broad spectral width of o and p resonances, we have 
folded the pion thermal width T^^m) by the vacuum spectral functions of those 
resonances as we have done in our previous work 17 to calculate shear viscosity. The 
multi-peak structure of total Y n u M ^{Mk) (solid line in the lower panel) is polished 
to a smooth curve (solid line in the upper panel) after introducing the folding effect. 
This is because the region of branch cuts are generally overlapped to each other 
during the folding operation, which is reflected by the dotted and dashed lines for 
7T(7 and np loops in the upper panel of Fig. 0- The on-shell contribution of the 
pion thermal width (i.e. r v (^M){Mk = m„)) has been indicated by the straight 
dotted line in the Fig. 0, where we see that the contribution of np loop dominates 
over the no loop and its total value becomes little smaller due to folding effect. 

This on-shell contribution from the Landau cut is associated with the forward 
and reverse scattering processes.^! In the forward process, the propagating n + may 
disappear by absorbing a thermalized n~ from the medium to create a thermal- 
ized p° or o. Whereas in the reverse process, the n + may appear by absorbing a 
thermalized p° or o from the medium as well as by emitting a thermalized n~. 

In upper panel of Fig. 0, the on-shell contribution of 7rcr (dotted line), np 
(dashed line) loops and total (solid line) thermal width for pion are plotted against 
temperature axis. Contributions of both loops are monotonically increasing func¬ 
tions but np takes leading part in the total contribution. The corresponding results 
for mean free path, defined by A n ( n M) = k/(u)^T n ^ wM ))i are shown in the lower 
panel of Fig. (fUj). Being inverse of thermal width, the mean free paths become 
decreasing functions of T. 

The momentum distribution for different components of thermal widths and 
mean free paths are presented in the upper and lower panels of Fig. 0. Upper 
panel shows that the pion thermal width for no (dotted line) and np (dashed 
line) loops become significant in the low and high momentum regions respectively. 
The corresponding mean free paths in the lower panel of the figure expose this 
complementary feature of two resonances in more prominent way. Here we notice 
that mean free path contribution from no loop is tending to diverge after k ~ 0.1 
GeV while divergent nature for np loop contribution is starting below the k « 0.2 
GeV. A phenomenological lesson from this fact is that the low momentum pion will 
get the relevant dissipation by scattering with o resonance whereas the dissipation 
of high momentum pion will be coming from the np scattering in the medium. 

Using the thermal width Y n ^M)(k,T) in Eq. 0, one can estimate the ther¬ 
mal conductivity, which possesses the temperature dependence from not only the 
thermal width but also the Bose-enhanced phase space factor of pions. The contri¬ 
butions of no (dotted line), np (dashed line) loops and their total (solid line) in 
the thermal conductivity are presented as functions of temperature in the Fig. 0, 
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Fig. 5. The contributions of ncr (dotted line), np (dashed line) loops and their total in the thermal 
conductivity k, are plotted against T axis. The upper and lower panels of this figure and latter 
Figs. m and 0 are allotted for demonstrating with and without folding results respectively. 



T (GeV) 

Fig. 6. T dependence of n without (dashed line) and with (solid line) form factor. 


where the upper and lower panels show the results with and without folding re- 
spectively. The k w (T) due to np loop appears as a decreasing function and diverges 
in low temperature region, where 7rer loop has a finite contributions. Because of 
this compensating nature of two loops, n^(T) for the mesonic system is revealed 
as non-divergent well behaved function in the entire temperature domain. Though 
the total thermal conductivity without folding effect (lower panel) is exhibiting a 
divergent nature at very low temperatures (T < 0.020 GeV) but this divergence is 
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Fig. 7. Total k(T) for different set of mass m. n (GeV) and width Y n (GeV) of a meson, taken from 
RefsHmmmSby following the previous Ref023 


cured completely after incorporating the folding effect (upper panel). 

Similar to folding effect, our results can also be influenced by cut-off parameter A 
of the hadronic form factor, which is taken in the ttkM vertex. The results without 
form factor (dashed line) yield the minimum values of re vs T curve, whereas results 
with form factor (solid line) exhibit the enhanced values of n as shown in Fig. 0. 
Owing to this fact, the thermal conductivity in our effective hadronic model may 
vary within a band; e.g. k = 0.0086 — 0.013 GeV 2 at T = 0.15 GeV. 

Besides this, our results may suffer from phenomenological uncertainty for fixing 
the parameters of a resonance as discussed in our earlier work on shear viscosity? 1 '’ 
Longstanding controversies about the properties of a resonance seem to be settling 
to a consensus 2 ® as the band of its mass-width values have been shortened from 
(ttv = 0.4 - 1.2 GeV, T a =0.6-1 GeV)P to (m CT = 0.4 - 0.55 GeV, T * = 
0.4 — 0.7 GeV)P^ Using the minimum (dash-dotted line) and maximum (dotted) 
values of (m a ,T a ) from latest PDGp^ the band of k without (lower panel) and 
with (upper panel) folding effect are shown in Fig. [3 The results of k, for the 
mass-width parameters taken from BEfpSED (solid line) and E791^ (dashed line) 
experiments, are also included in the Fig. 0. Among them, the set of (m CT , IV) from 
BES experiment is arbitrarily chosen to generate all the previous results, shown in 
Fig. 0 to j6|). The list of effective coupling (g a ) is provided in the Table of RefP^ 
From the upper panel of Fig. 0, n at T = 0.15 GeV may vary from 0.012 GeV 2 to 
0.016 GeV 2 . Hence, adding the uncertainty of a meson parameters with the form 
factor effect, we get n(T = 0.15 GeV) = 0.0086 — 0.16 GeV 2 , which is exhibiting 
quite large band of prediction of thermal conductivity from the effective hadronic 
model. From the Table 0, one can notice that the numerical values of k of earlier 
RefsPHEE] are within the band of uncertainty, contained in our results. 
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4. Summary 

In Summary, an estimation of thermal conductivity for pionic medium has been 
done with the help of the effective hadronic model. We have started with the stan¬ 
dard expression of thermal conductivity in RTA approach, where inverse of pion re¬ 
laxation time (r n ) or pion thermal width r„. has been evaluated from the imaginary 
part of pion self-energy. Keeping in mind about the appearance of a, p resonances 
in the ttw scattering cross section, we have evaluated pion self-energy for na and 7r p 
loops in the real-time formalism of thermal field theory. We get a non-monotonic 
momentum distribution of pion thermal width, which has to be integrated out by 
the Bose-enhanced phase factor during estimation of thermal conductivity. Our es¬ 
timation has a large band of uncertainty because of the various phenomenological 
parameters such as mass-width parameter of sigma meson, hadronic form factor etc. 
SomePsMni of the earlier estimation of thermal conductivity for strongly interact¬ 
ing matter in hadronic temperature domain are within our band of estimation. 
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